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Abstract 



Q^i In this paper, we report an interesting integrable equation that has both solitons and 

^ ■ kink solutions. The integrable equation we study is (— ^|^)t = 2uux, which actually 

comes from the negative KdV hierarchy and could be transformed to the Camassa- 
Holm equation through a gauge transform. The Lax pair of the equation is derived 
to guarantee its integrability, and furthermore the equation is shown to have classical 
' solitons, periodic soliton and kink solutions. 

o 

\0 ■ PACS. 02.30.Ik C Integrable systems. 

PACS. 05.45.Yv C Solitons. 

PACS. 03.75.Lm C Tunneling, Josephson effect, Bose-Einstein condensates in periodic po- 
tentials, solitons, vortices, and topological excitations. 

> : 

H : 1 Introduction 

Soliton theory and integrable systems play an important role in the study of nonlinear water 
wave equations. They have many significant applications in fluid mechanics, nonlinear optics, 
classical and quantum flelds theories etc. Particularly in recent years, more focuses have been 
pulled to integrable systems with non-smooth solitons, such as peakons, cuspons, since the 
study of the remarkable Camassa-Holm (CH) equation with peakon solutions [2]. Henceforth, 
much progress have been made in the study of non-smooth solitons for integrable equations 
[4-20, 22, 27-33]. 

In this paper, we consider the following integrable equation 

{^)t = 2uux, (1) 
u 
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which is actually the first member in the negative KdV hierarchy [29] . Equation ([T]) is proven 
equivalent to the Camassa-Holm (CH) equation: 'mt + mxU + 2mUx = 0,m = u — u^x through 
a gauge transform (see Remark 1 in the paper). Therefore, we find a simpler reduced form of 
the CH equation. The Lax pair of the equation ([1]) is derived to guarantee its integrability, 
and furthermore the equation is shown to have classical solitons, periodic solitons and kink 
solutions. 

2 Derivation of Equation (jlj) and Lax Representation 

Let us consider the Schrodinger-KdV spectral problem: 

Li) = i^xx + = \i), (2) 

where A is an eigenvalue, ip is the eigenfunction corresponding to the eigenvalue A, and v is 
a potential function. One can easily get the following Lenard operator relation: 

KVX = AJVA, (3) 

where VA = |^ = ?/;^ is the functional gradient of the spectral problem ([2]) with respect to v, 
K = + \ {vd + dv) and J = d are two Hamiltonian operators as known in the literature 

By setting v = — we have the product form of operators K, L and their inverses 

K = lu-^du^du^du-^, K-^ = iu^d'^u-^d'^u^^d'^u^, 
L = ^d~^u~'^dv?dv?du~'^ , L^^ = Av?d~^u~'^d~^u~'^d~^v?d^ 
L = d'^ + V = u^^du^du"^ , L^^ = ud^^u^'^d^^u, 

where £ = J~^K and its inverse £~ = K"^J are the recursion operators for the positive 
order and negative order KdV hierarchy that we study below. 

Now, according to the Lenard's operators K and J, we construct the entire KdV hier- 
archy, and then we show the integrability of the hierarchy through solving a key operator 
equation. 

Let Go e Ker J = {G e | JG = 0} and G_i e Ker K = {G e G°°(M) | KG = 0}. 

We define the Lenard's sequence 

(4) 

where L,L~^ are defined by Eq. ([2]). Therefore we generate a hierarchy of nonlinear evolution 
equations (NLEEs): 

vt, = JGk = KGk-i, yk G Z, (5) 
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which is called the entire KdV hierarchy. We will see below that the positive order [k > 0) 
gives the regular KdV hierarchy usually mentioned in the literature [T], while the negative 
order [k < 0) produces some interesting equations gauge-equivalent to the Camassa-Holm 
equation [2]. Apparently, this hierarchy possesses the bi-Hamiltonian structure because of 
the Hamiltonian properties of K, J. Let us now give special equations in the entire KdV 
hierarchy (jS]). 

• Choosing Gq = 2 E Ker J (therefore Gi = u) leads to the second positive member of 
the hierarchy (E]): 

13 

Vt2 = -j^xxx + -^'W^, (6) 

which is exactly the well-known KdV equation. Here there is nothing new. Therefore, 
the positive order (/c > 0) in the hierarchy (jS]) yields the regular KdV hierarchy usually 
studied in the literature [1]. 

• Now, let us find kernel elements G„i G Ker K in order to get the negative member 
of the hierarchy ([5]). Due to the product form of K and K~^, G^i = K~^0 has the 
following three seed solutions: 

Gi, = f{tnW, G\ = g{tnWd-^u-\ G^_^ = h{tnWd-^u-^d-^u-\ 

where f{tn),g{tn)ih{tn) are three arbitrarily given functions with respect to the time 
variables tn, but independent of x. They produce three iso-spectral (Aj^. = 0) negative 
order KdV hierarchies of Eq. (j5]) 

vt, = JL'^' ■ Gi,, / = 1, 2, 3, k = -l, -2, .... (7) 

When k = —1, their representative equations are: 

(-— ) = 2f{Quu^, (8) 

'-) = g (tn) {2uu^d~\-^ + 1) , (9) 



f-^) = h{t^)i2uu^d-^u'^d~^u~^ + d-^u~^). (10) 

Remark 1. Apparently, the first one is differential and simpler and exactly recovers the 
equation ([1]) after setting f (fn) = 1, which we focus on in the current paper. Actually, these 
three representative equations and / flOj) come from Vt_^ = JG-i = JK~^0 with 

V = — Clearly, Vt_^ = JK~^0 is equivalent to KJ'^Vt = (t = t_i), that is, 

'^txx 



/ x V t>^ / : 



2vt = 0. (11) 



This equation is exactly the one studied by Fuchssteiner ITSf (see equations (7.1) and (7. 
there. Equations (7.1) in ITSf has a typo and should be same as 01]) . From ITS^ . the 
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Camassa-Holm (CH) equation is gauge- equivalent to equation [Jl]) through some hodograph 
transformations (7.11) and (7.12) in JWI. In our paper, through using v = we further 

reduce equation l[TJ\) to a more simple form (i.e. equation i^^): 

^""^ 2uu^. (12) 



u / t 



In other words, we found a very interesting fact that equation (fT2l) can be viewed as a reduc- 
tion form of the CH equation due to the above gauge- equivalence. In next section, we will 
solve this form. 

In the paper [19], the author deah with equation (9^ + 4v + 2vxd~^)vt = by using 
the positive KdV hierarchy approach, and all soliton solutions were given implicitly. This 
equation could be transformed to {—Uxx/u)t = 2uu.j. through v = —u^x/u, like we mentioned 
earlier in our paper, but, this is only one of three reductions. So, solutions of this equation 
can not give all solitons of our equation {—Uxx/u)t = 2uUx. In our paper, we present all 
solitons and kink solutions in a explicit form. Also, there is the connection of the first 
negative KdV equation (9^ + Av + 2vxd~^)vt = with sine-Gordon [20]. But, the equation 
{—Uxx/u)t = 2uUx we propose in the current paper is not equivalent to the sine-Gordon 
equation, because the sine-Gordon equation has only kink solution while our equation has 
both kink solutions and classical solitons. 

Of course, we may generate higher order nonlinear equations by selecting different mem- 
bers in the hierarchy. In the following, we will see that all equations in the KdV hierarchy 
([5]) are integrable. Particularly, the above three equations ([8]), ([9]), and (ITOl) are integrable. 

Let us return to the spectral problem ([2]). Apparently, the Gateaux derivative matrix 
L*(,^) of the spectral operator L in the direction ^ G C°°(M) at point v is 



f/(n + eO=e (13) 



e=0 



which is obviously an injective homomorphism, i.e. t/*(^) = <(=^ ^ = 0. 

For any given C°°-function G, one may consider the following operator equation |3] with 
respect ioV = V{G) 

[V, L] = (KG) - U{JG)L. (14) 

Theorem 1. For the spectral problem ([2]) and an arbitrary -function G, the operator 
equation has the following solution 

V = -\Gx + \Gd, (15) 
where d = dx = -§;^, and subscripts stand for the partial derivatives in x. 
Proof: A direct substitution will complete the proof. 
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Theorem 2. Let Gq E Ker J, G^i G Ker K, and let each Gj be given through the Lenard 
sequence Then, 

1. each new vector field = JG^, k E Z satisfies the following commutator representa- 
tion 

L,{Xk) = [Vk,L], \/keZ; (16) 

2. the entire KdV hierarchy Ij^, i.e. 

vt, =Xk = JGu, yk e Z, (17) 
possesses the Lax representation 

Lt^ = [Vk,L], VfcGZ, (18) 

where 

E k > 0, 

V, = J2v{G,)L^'-^-'\ 5^ = <! 0, k = 0, (19) 



and V{Gj) is given by Eq. ( f75j) with G = Gj. 

Proof: Let us only prove the case for A; < 0. We have 



-1 



[Vk,L] = -J2[V{G,).L]L'^~^~' 

j=k 
-1 

= - ^(L,(A'G,) - L,{JG,)L)L''-^-^ 



j=k 

= -J2MKGj)L'^''^'^^ - L,{KGj.i)L''^' 

j=k 

= L,{KGk^,) - L,{KG.,)L' 
= L,{KGk^i) = L,{JGu) 
= L^{Xk). 

Noticing = L^{vtJ, we have 

Lt,-[Vk,L] = L.ivt,-X,). 

The injectiveness of L^, implies the second result holds. 

So, the entire KdV hierarchy has Lax pair and all equations in the hierarchy are 
therefore integrable. In particular, the KdV equation (|6]) has the Lax pair Lt^ = [Wi,L] 
with L = d"^ + V and Wi = + ^vd + jVx, which was well-known in the literature [I]. 
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Interesting thing is that the first member [k = —1) in the negative order KdV hierarchy 
([7]) has the standard Lax representation L^ ^ = [V"1]^,L] with V"i^ = (^CLi^; — ^GLiS) L~^, 
I = 1,2,3, L = d"^ + V = u~^dv?du~^, and L"^ = ud~^u^^d~^u. In particular, the negative 
KdV equation ([8]) possesses the following Lax form: Lt ^ = with 

= ( -uu^ u^d^L~^ 

\2 2 J 

= — uo u. 
2 

All of those negative members in the hierarchies ([7]) are integrable. 



3 All traveling wave solutions of ([I 



Let us now consider the traveling wave solution of equation ([T]) through a generic setting 
u{x,t) = U{x — ct), where c is the wave speed. Let ^ = x — ct, then u{x,t) = U{^). 
Substituting it into equation ([T]) yields 



c{'L)' = 2UU'. 



(20) 



Integrating it once, we obtain the following standard cubic Hamiltonian system for c 7^ : 

dH 



dy'd^ c 



dU' 



where g is an integral constant, and the Hamiltonian function is 



(21) 



(22) 



When gc > 0, (21) has only one equilibrium point 0(0,0). When gc < 0, (21) has three 



equilibrium points 0(0,0) and Ei^2i^\'\cg\,0). Write that 

ho = H{0, 0) = 0, hi = H{±^/\^\, 0) = ^cg''. 

By qualitative analysis, we have the bifurcations of phase portraits of (21) in the 
(c, (7)— parametric plan shown in Fig.l (l-l)-(l-6). 




(1-1) 



;i-2) 
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;i-3) 



(1-4) (1-5) (1-6) 

Fig. 1 The change of phase portraits of (21) in the (c, (?)— parameter plan. 

(1-1) .g = 0,c > 0. (1-2) ,9 > 0,c > 0. (1-3) g > 0, c < 0. (1-4) .g = 0, c < 0. (1-5) .g < 0,c < 0. (1-6) 
g < 0, c> 0. 

Next, we present the exact travehng wave solutions of (1) in an explicit form. 
Case 1: c>0,g = (see Fig.l (1-1)). 

In this case, corresponding to the saddle point 0(0, 0) (12T!) reads as y = Using the 

first equation of (21) and taking integration, we obtain 



^ = X — ct, 



(23) 



where .^o is an initial value of ^. Clearly, when ^ — )■ — .^o, U{S,) — t- oo. i.e., U{S,) is unbounded 
at ^ = ^0- Thus, we have two unbounded breaking wave solutions shown in Fig. 2. 




Fig. 2 The profiles of the functions (23) where = ■pl- 



ease 2: c> 0,^ > (see Fig.l (1-2)). 

Corresponding to the saddle 0(0,0) ( 1211) reads as = 

gU^ + yU*- Using the first 



equation of (21) to take integration, we obtain 



U{0 



8Acg 



(24) 



where A is an integrant constant. When A"^ = 8cg, the functions defined by (24) are 
discontinuous at ,^ = — ,^o- The profiles of (24) like Fig. 2. 
Case 3: ^ > 0,c < (see Fig.l (1-3)). 

There exist three equihbrium points of (21) at £'1,2 and 0(0, 0). O is a saddle point, £'1,2 
are center points. 
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Corresponding to two homoclinic orbits, defined by H{U, y) = 0, we have the parametric 
representations 

f/(0 = iV^sech^MZ^. (25) 

They give two sohton solutions of (1). 

For h E {jcg'^,0), corresponding to two families of periodic orbits of (21), defined by 
H{U,y) = h, i.e., y ^ = ^(4|c|/i + 2g\c\U' - U^) = - U^){U^ - rl), where r? = 

g\c\ + \/ g^c? — ^\c\h^r\ = g\c\ — \/ g^c? — 4|c|/;,, we obtain the parametric representations of 
periodic wave solutions of (1) as follows: 



^(0 = ±ridn I ^^^^^ ) . (26) 

For h G (0, oo), corresponding to the family of periodic orbits of (21), enclosing three 
equilibrium points defined by (f/, y) = h, we have the following parametric representation 
of periodic wave solutions of (1): 



-(e) = ricnh/-^e,^^^l. (27) 



Case 4: ^ < 0,c < (see Fig.l (1-4), (1-5)). 

In this case, the origin 0(0,0) of (21) is an unique equilibrium point, which is a center. 
There exists a family of periodic orbits of (21), enclosing the origin. (1) has the same 
parametric representation of periodic wave solutions as (27). 

Case 5: ^ < 0, c> (see Fig.l (1-6)). 

In this case, there exist three equilibrium points of (21) at £"1,2 and 0(0, 0). O is a center, 
£^1,2 are saddle points. The heteroclinic orbits, defined by H{u, y) = hi, have the parametric 
representations 

f/(0 = ±v^tanh (^^^ j > (28) 

which gives a kink wave solution and an anti-kink wave solution of (1). 

We see from (22) that y'^ = ^{4:ch + 2cgU'^ + Lf^. For h G (0, hi), it can be written as 
= Yc^i^l ~ ^2)(2| - f/^), where zl = \g\c+ g^c^ - Ach, z\ = \g\c- ^/g'^c^ - 4:ch, Thus, 
the family of periodic orbits, defined by H{u, y) = h, has the parametric representation 

u{0 = Z2sn ^] , (29) 



which gives a family of periodic wave solutions of (1). 

4 Conclusions 

In this paper, we reported an interesting property of integrable system: solitons and kink 
solutions can occur in the same integrable equation, and those solutions are given explicitly. 
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Within our knowledge, this is probably the first integrable example possessing such property. 
We found this equation in the negative order KdV hierarchy, which is gauge-equivalent to 
the CH equation. Since equation ([T]) has the Lax pair, we may try to get r-matrix structure 
of the constrained system of Lax equations, and parametric and algebro-geometric solutions 
[28] . but that is beyond the scope of this paper. The symmetry of equations (8), (9), and 
(10) were already discussed in [23]. Recently, a two-fold integrable hierarchy associated with 
the KdV equation was given in [21]. About other negative order integrable hierarchies, such 
as the AKNS, the Kaup-Newell, the Harry-Dym , the Toda etc, one may see the literature 
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